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We recently determined the mass of the most massive star known to the date, R136al with a mass 
at birth 320 times the mass of our sun, as well as the mass of several other stars that are more 
massive than 150 Mq. Such massive stars (~ 150-300 Mq) may end their life as pair-instability 
supemovae (PISN) if they retain enough mass until they die. We have calculated a grid of stellar 
evolution models in order to investigate the impact of mass loss and rotation on the evolution 
and fate of these very massive stars. As mass loss is very strong at solar metallicity, our models 
predict that most of the very massive stars will die as type Ic SNe. Only slowly and non-rotating 
stars at metallicities below that of the LMC might retain enough mass to produce a PISN. This 
would mean that the first stellar generations might have produced PISN although their chemical 
signature is not observed in extremely metal poor stars in the halo of our galaxy. 
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1. Introduction 

The fate of a star is determined by its mass, composition, rotation rate and binarity, the mass 
being the most important parameter. Massive stars with an initial mass between 8 and 140 
will evolve through all burning stages and form an iron core. The fate of the star follows this 
sequence as a function of increasing initial mass: SNII-SNIb-SNIc-BH without SN. For initial 
masses between 140 and 260 M©, the stars are thought to die as pair-instability supernovae if they 
do not lose much mass during their evolution, which is expected for low-metallicity non-rotating 
stars [1, 2]. PISNe are very bright and produce large quantities of iron [1]. Even if they are 
rare, they contribute significantly to the chemical evolution of galaxies. Their existence is debated 
because their chemical signature is not observed in extremely metal poor stars [3]. The absence of 
this chemical signature could be due to the fact that stars more massive than 140 do not form. 
This argument was supported by local observations. Indeed, the initial mass function (IMF) of the 
most massive young cluster in the Galaxy measured using HST seems to indicate a lack of stars 
with mass more massive than 150 M© [4]. However, observations of 2007bi [5] indicate that it is 
consistent with PISN theoretical models. 

In this work, we present evidence supporting the existence of stars more massive than 150 
M0. We also investigate the effect of rotation, mass loss and metallicity on the evolution of very 
massive stars. Main sequence evolutionary models are compared with the observations and we find 
an excellent agreement for the NGC 3603 and R136 systems [6]. Finally, we discuss the fate of 
these very massive stars. 

2. Stellar Evolution Models 

A grid of stellar evolution model for 85, 120, 150, 200, 300 and 500 M© stars has been 
calculated until at least the end of the main sequence using the Geneva stellar evolution code [7] 
both with and without rotation. These models include the physics of rotation and mass loss, which 
are both very important in the evolution of very massive stars. The models include all the effects of 
the rotation: centrifugal support, mass-loss enhancement and especially mixing in radiative zones 
[7]. In these models, we compute the main sequence evolution using theoretical mass loss rates 
for O type star [8], which match observations very well [9]. We consider that the star enters the 
Wolf-Rayet phase when the surface hydrogen content Xh < 30% if Teff > 10,000 K, during which 
an empirical mass-loss calibration is followed [10]. We choose the ratio of the initial velocity to 
critical rotation of Vinu/vcrit^ 0.4 for the rotating models, which corresponds to surface equatorial 
velocities around 350 km~^ for 85 M© model and 450 km~^ for the 500 M© model. 

The evolution tracks of 120 M© models are presented in Fig. 1. In the figure we can see 
that the effect of rotation is significant. Additional mixing above the H-burning convective core 
causes helium to be mixed up to the surface. Thus the opacity in the outer layers decreases and 
the surface of rotating models remains hotter than that of non-rotating models. During the main 
sequence, the effective temperature of the rotating models remains hot, around 45 000-55 000 K 
whereas for the non-rotating models the effective temperature decreases steadily to 20 000-25 000 
K. Hence, rapidly rotating stars evolve directly to the classically W-R phase while the slow rotator 
stars are expected to become T] Carinae-like blue luminous variable stars (see Meynet and Meader 
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Models: l:H120z06S400, 2:H120z06S000, 3:H120zl4S000 
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Figure 1: Comparison between main sequence evolutionary predictions for both rotating (Z = 0.006: 
H120z06S400) and non-rotating (Z = 0.006:H120z06S000, Z = 0.014: H120zl4S000) 120 M© models. 
For model H120z06S400, the horizontal lines in the top left-hand panel correspond to the transition to the 
WR phase, phase during which the photosphere is in the wind rather than at the surface of the star. 



2005 [11]). In Fig. 1, we compare solar metallicity with LMC metallicity models and we can see 
that lower metallicity models reach higher luminosities. This is due to the lower mass loss. 



3. Mass Determination of the Most Massive Stars 

The results from the evolutionary models were compared with spectroscopic analysis of two 
young star clusters, i.e. NGC 3603 and R136. Details of the spectroscopic analysis of the star 
clusters can be found in Crowther et al. [6]. In Fig. 2, we compare the evolutionary models of 
solar metallicity with the observational properties of NGC 3603 WN stars while in Fig. 3, the 
evolutionary models of LMC metallicity are compared to the properties of R136 stars. 

For NGC 3603, we find a good agreement with dynamical masses for Ala and Alb for ini- 
tially non-rotating models at ages of 1.5 ± 0.1 Myr. This provides strong support for our mass 
determination method (for more details see Crowther et al. 2010 [6]). For R136, we obtain good 
agreement between the R136 stars and initially rapidly rotating, 165 - 320 M© models at ages of 
^1.7 zh 0.2 Myr. In particular, we determine a birth mass of 320 M© for R136al, which is the most 
massive star known to date and the mass of which exceeds the commonly accepted upper limit 
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Figure 2: Comparison between Solar metallicity (Z = 1.4%) models calculated for the main-sequence evo- 
lution of 85 - 200 M0 stars (initially rotating at Vmit/vcrit =0-4 [dotted] and [solid]), and the physical 
properties derived from spectroscopic analysis of NGC 3603 WN6h stars. 



of I]V[F significantly. Note also that observed mass-loss rates match to within 0.2 dex theoretical 
predictions [8] at solar metallicity and LIMC metallicity for NGC 3603 and R136, respectively. 



4. Fate of Very IMassive Stars 

The fast rotating very massive stars around solar metallicity experience strong mass loss and 
are thus predicted to end their lives as type Ic SNe. IVleanwhile for the slow rotating and non- 
rotating stars, the mass loss effect is less than for the rotating stars. In Table 1, we present the 
hydrogen surface content for 120, 150 and 200 Mq non-rotating stars at different metallicities at 
the end of the main sequence. Although the models will be evolved further in order to obtain a 
more precise picture of the fate of these stars, if the star has lost all of its hydrogen envelope by 
the end of its main sequence phase, it will become a WR star and probably die as as type Ic SN. 
This is the case for the solar metallicity models. Stars at SIMC metallicity, on the other hand, retain 
some of their hydrogen rich envelope at the end of the main sequence. These stars will go through 
the LBV phase, where the mass loss is poorly understood. They may still retain enough mass to 
die as PISN. For stars at LIMC metallicity the fate depends strongly on the mass loss and the initial 
mass of the star. The fate of very massive stars at lower metallicities will be investigated in a future 
work. 
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Figure 3: Comparison between LMC-metallicity models calculated for the main-sequence evolution of 85 
- 500 M0 stars, initially rotating at vinit/vcrit = 0-4 (dotted) or (solid) and the physical properties derived 
from our spectroscopic analysis. 
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Table 1: Surface of hydrogen content predicted at end of the main sequence for the 120, 150 and 200 
non-rotating stars. 
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